Effect of diffusion of elements on network topology and self-organization 
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ABSTRACT 

We study the influence of elements diffusing in and out of a network to the topological changes of 
the network and characterize it by investigating the behavior of probability of degree distribution 
(F(fc)) with degree k. The local memory of the incoming element and its interaction with the 
elements already present in the network during the growing process significantly affect the network 
stability which in turn reorganize the network properties. We found that the properties of T(k) 
of this network are deviated from scale free type, where the power law behavior contains a 
exponentially decay factor supporting earlier reported results of Amaral et.al. [ll[ and Newman 
12| and recent statistical analysis results on degree distribution data of some scale free network 
131 ] - Our numerical results also support the behavior of this T(k). However, we found numer- 
ically the contribution from exponential factor to the r(fe) to be very weak as compared to the 
scale free factor showing that the network as a whole carries the scale free properties approximately. 

KEYWORDS : Scale free network, Probability of degree distribution, Master equation, Preferential 
attachment, Growing network. 
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Complex networks in which elements interact among 
each other is recently focused on because these networks 
fall on mostly in real and natural networks such as scale 
free and small world networks 0-H[,_both are microscopic 
and macroscopic natural networks [HEIQ- Further these 
networks possess growth and high clustering nodes, which 
are some of the most important properties of real natural 
networks, due to activation and deactivation of the nodes 
in the network 0, @, • If new element is introduced in 
the network, it may be added or rejected that leads to 
local perturbation in the network topology which in turn 
affects in the network stability locally as well as globally 
• This effect of local stability on the network topology 
gives rise correlation between network stability and the 
emergence of self-organization within the network itself 

0- 

The investigation of properties such as self- 
organization, network stability, signal transduction 
etc in biological networks in particular has been im- 
portant area of research in recent times 0. In these 
networks, where biochemical reactions are considered 
to be links among the various participating molecules 
which are taken as nodes, diffusing in and out of intracel- 
lular (various biomolecules) and intercellular molecules 
(signaling molecules) via intra and extracellular media is 
an important common phenomena . As a consequence 
of the diffusing molecules, activation and deactivation 
of nodes in the network take place that lead to change 
in the network topology and information transduction 
within the system and among the systems, where scale 
free and small world effects ( pre sence of highly clustering 
nodes) are the signatures [{J, [l(| • 

However, it has been pointed out that many of the real 
networks do not follow scale free nature since the degree 
distributions of these networks do not follow power laws 



but associate with certain exponential component in the 
distribution function [ll|, E3 • ^ * s a l so argued that if one 
does proper statistical analysis of the degree distribution 
data instead of simply fitting the data by straight line in 
the log- log plot, the degree distribution function follows 
power law within a certain small range defined by a cut 
off parameter, k c and beyond this value it follows a sharp 
drop-off given by a exponential factor in the distribution 

El. 

The growth dynamics of a complex network can be de- 
scribed by the change in directed links among the nodes 
due to the fluctuation in the network topol ogy by incom- 
ing and outgoing elements in the network [12| ■ Consider 
ki be the number of links incident on a particular node 
"i" or indegree of node "i". Because of this incident 
node, there should be local perturbation in and around 
node "i" which in turn affect the network topology. So 
there could be two possible states of each node : active 
and inactive states (l2| . The added new node is consid- 
ered to be always in the active state first and recieves 
new links from the generated node until it gets deacti- 
vated 0, G3 ■ Once a particular node is deactivated, it 
will no longer recieve links anymore. Now consider P(k) 
as the probability of deactivation rate which decreases 
as in-dcgree node increases. If we assume this P to be, 
P{k) ~ constant + f (k) , its specific form is considered to 
be given by the following form, 

1 ib 

P{k)-- + ^-T (1) 
e 77 + k 

where, e, tp and 77 are constants provided e, 77)0. The first 
term in the equation ([1} is a constant (^-(1) which indi- 
cates that the deactivation rate was initially at constant 
in the network. The second term incorporates the pertu- 
abation in the deactivation probability when new nodes 



are added. 



A. Estimation of probability of degree distribution 

Now if we take F(fc, t) to be the in-dcgree distribution 
of the active nodes at any instant of time t, then the time 
evolution of T(k, t) is determined by the following Master 
equation (i~2"j . 



r(fc + i,t + 1) = [i -p(k)}r{k,t) 



(2) 



Using the equation ((T|) in $2$, the following equation can 
be obtained, 



T(k+l,t + l)-T(k,t) 



>P 



V 



T(k,t) (3) 



Applying boundary condition, T(0) = constant which 
means that constant rate of new nodes. The fluctuation 
is assumed to be small and taken to be constant. Now 
keeping time to be approximately constant, equation ([3]) 
can be written as, 



r(fc + i)-r(*) 



1 ip 
- + — - 

e i] + 



T(k) 



(4) 



Now, taking k to be continuous, Master equation ((4]) is 
reduced to the following differential form, 



dk 



V 



T(k,t) 



(5) 



Solving equation ([5]), we obtained the expression for T(fc), 

r(jfe) = /3 x e-^in + k)-^ (6) 

where, j3 is a constant. Since the total number of nodes 
is large as compared to the number of active nodes, the 
overall degree distribution, W(k) can be approximated 
by considering active nodes only, which is given by, 



W{k) = 



dT(k) 
dk 



ae 



(v + k) 



(7) 



where, a is a constant given by, a = fiipe^ and tpe))(n + 
k). The constant (3 in equation (O can be calculated by 
evaluating the normalization condition i.e. J Q T(k)dk = 
1, giving the value, /3 = (n + . Again solving the first 
moment of the active node distribution, 



kW{k)dk 



(8) 



we can estimate the exponent ip. Now defining ip = 7 — 1, 
where 7 is the constant related to normalization factor of 
the probability T(fc) which we obtain after simplification 

as, 



7 = 1 + q(m, rj, e) 



(9) 
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FIG. 1: Plot of probability of in-degree distribution of active 
nodes T(k) as a function of k : (a) T(k) vs k for different 
values of e; (b) T(fc) vs k for different values of ip. The values 
of the constants are: rj = 10, ft — 20. 



Hence, 7 depends on m, rj and e. 

Now, we plot r(fc) in equation (JB]) as a function of k 
for different values of e and ip as shown in upper and 
lower panel of Fig. 1. The plots indicate that for all 
values of e, the signature of exponential component re- 
flects only in the small values of k with small range. In 
the remaining large range of k, the T(fc) is dominated by 
power law behavior (straight line nature of lnT(k)-ln{k) 
data). Moreover, as the value of e increases, power law 
part dominate over exponential component in a wider 
range of k. 

Now we try to estimate the order parameter ip from 
the plot of r(fc) as a function of k for different values of 
i/j shown in Fig.l lower panel. The plot shows that ip is 
strong for scale free behavior when ip < 2.5. These char- 
acters of ip and e indicate that the network of this type 
has weak contribution in T(fc) from the exponential part, 
but scale free factor dominates wider range of k. This 
indicates that the overall network properties are domi- 
nated and controlled by scale free properties supporting 
Perotti et al 0). 
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FIG. 2: Plot of InT(k) of active nodes as a function of ln(k) for three different values of 6 = 2, 4, 6, N = 200, 500, 1000 and 
fixed value of q — 3. The curve lines are fitted curves on the data by InT(k') = A — Bk' — Ce Dk . The dotted straight lines are 
the fitted lines by straight line IriT(k') = A — Bk' to the certain range of data points. 
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N N 

FIG. 3: Plots of different parameters as a function of N : (a) 
i\) versus N, (b) D versus N, (c) C versus N and (d) %j) s j 
versus N. The error bars are coming from the calculation of 
different parameters from various data sets i.e. various data 
of b, q and N. 



B. Preferential attachment 

If A(k) is the degree dependent attachment rate and 
the network contains t nodes, then tW(k) of these will 
have degree k. The number of active nodes with degree 
k is ?7ir(fc). Hence, the average increase in degree is 
given by, A(fc) = . Using equation ©, we get the 

preference attachment rate as, 

A(/c) = Z(r) + k)e~^ (10) 
where, the constant Z is given by, Z = j£. 

C. Numerical results 

We follow the numerical algorithm due to Pcrotti ct al 
to calculate probability of degree distribution, T(k) of 
growing network generated by diffusion of nodes in and 
out of the network. The algorithm is based on generated 
memory of grow or shrink at every step characterized by 
eigen value (A) of the modified Jacobian matrix of the 
elements in the network at that step. If A ( then the 
new node diffused inside the network is accepted by intro- 
ducing a range of interaction range defined by a uniform 
random number b in [—6, b]. However the condition A ) 
destablize the network, a condition is imposed that ei- 
ther the node is eliminated or certain number of already 
existed nodes (q < in the previous step are elimi- 

nated. 

Now we present our numerical results based on the 
algorithm we explained for different values of b = 
1, 2, 4, 6, 8, 10 and q = 2, 3, 4 as shown in Fig.2. We cal- 
culated r(fc) as a function of k for the mentioned param- 



eters and logT(k) — log(k) for b = 2,4,6 for q = 3 is 
shown in the plots. The the curve lines in each plot are 
the fitted curves using our simplified expression of T(k) 
in equation ^ (we take 77 = for the sake of simplic- 
ity)!^. InT(k') = A-Bk'-Ce Dk ' , where, k' = ln(k) and 
A, B, C and D are constants to be obtained by fitting 
the curves. The equation fits well the data but the order 
parameter ip is found to be very small as compared to 
the predicted values of scale free networks. The param- 
eter ip is calculated function of network size (N), 
by averaging its values calculated for various values of 
parameters b and q, and is shown in Fig. 3 (left upper 
panel). Taking mean of the values for different values of 
b, q and N, we found V = 0.83 ± 0.053. Then we fit the 
curves with straight lines so that we can get the values 
of scale free parameter ip s f for the curves (fitted straight 
lines are shown in Fig.2) as a function of N. The result is 
shown in Fig. 3 (right lower panel) and the average value 
of ip sf is found to be 2.78 ± 0.127. 

We then search for the range of k within which scale 
free properties work by looking at the values of exponen- 
tial factor in the equation ©. So we obtain the values 
of C from the curve fitting procedure in all cases and is 
found to be very small and decreasing as N increases 
and become approximately stationary after N ~ 800 
(C = 0.47 x 10~ 6 ). The plot of C versus N is shown 
in lower left panel of Fig. 3. The scale free works in a 
small regime of k for all network sizes N, for example 
for N = 1000, it works within the regime 412 ( k ( 894 
approximately. 

We now search for the values of D as a function of N 
following the fitting process of the same set of data and 
is shown in Fig. 3 upper right panel. The average value 
of D is found to be 2.6 ± 0.256. Putting all these values 
of constants, we can write approximately the expression 
for T(fc) as in the following, 

r(jfe) RJ constant x e -4-7xi - 7 fc 2 6 x ^ + k yip ( U ) 

This means that the exponential factor influence the scale 
free properties slightly. The equation pip further shows 
that even if the scale free nature works in a small range of 
k, it dominates and carries the overall network properties. 
Eventually, the order parameter should scale by ip —> ip s f 
when one considers overall network topology. 

The activation and deactivation of nodes in a network 
due to diffusing of elements in and out of the network, 
provided the probability of deactivation rate we have 
taken in equation ([1}, definitely give a change in the 
in-degrec distribution of active nodes. This distribution 
contain a power in k factor with an exponential compo- 
nent that makes a different distribution from scale free 
type. When e))k, T(fc) follows power law falling to scale 
free behavior. But when e is comparable to k value, T(k) 
is deviated from scale free type. However the overall net- 
work carries the scale free properties as evident from our 
analytical results. 

Since the exponential factor contributes a little in the 
r(fc) as evident from our numerical results, the overall 
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network properties fall into scale free type. This indi- 
cates that the perturbation due to diffusion of nodes in 
and out of the network influnce the properties of net- 
work structure as indicated by the deviated value of or- 
der parameter ip. However still the scale free properties 
dominate the overall network properties indicating the 



self-organization of the network after the pertubation is 
induced. 
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